The Casimir interaction of a massive vector field between concentric spherical bodies 
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The Casimir interaction energy due to the vacuum fluctuations of a massive vector field between 
two perfectly conducting concentric spherical bodies is computed. The TE contribution to the 
Casimir interaction energy is a direct generalization of the massless case but the TM contribution 
is much more complicated. Each TM mode is a linear combination of a transverse mode which is 
the generalization of a TM mode in the massless case and a longitudinal mode that does not appear 
in the massless case. In contrast to the case of two parallel perfectly conducting plates, the are no 
TM discrete modes that vanish identically in the perfectly conducting spherical bodies. Numerical 
simulations show that the Casimir interaction force between the two bodies is always attractive. 
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I. INTRODUCTION 



Casimir effect has been under intensive investigation for its relation to many other areas of physics [1[ . It has been 
studied for both massless and massive scalar fields as well as spinor fields in different geometric configurations. For 
electromagnetic fields (massless vector fields), the situation is more subtle. Although Casimir effect of electromagnetic 
O j' fields has been extensively studied, there are very few works that considered Casimir effect of massive vector fields. 

The pioneering work of Barton and Dombey Q considered Casimir effect of a massive vector field on a pair of perfectly 
£^ ■ conducting plates in vacuum. This work revealed that the analysis for a massive vector field is much more complicated 
than for a massless vector field. One of the fundamental difference is that for a massless vector field, one can regard 
the electric field and magnetic field as the primary quantities and there is a gauge degree of freedom in the potentials. 
However, for a massive vector field, the primary quantities are the scalar and the vector potentials and the gauge 
freedom is lost. A massless vector field cannot penetrate through a perfectly conducting plate but for a massive vector 
field, although the electric field and the magnetic field have to vanish in a perfectly conducting plate, the potentials 
need not vanish. The latter implies that even for the Casimir effect of a massive vector field on a pair of perfect]; 
conducting plates, one has to treat it as if considering Casimir effect in dielectric plates. In fact, it was shown in 
that in the configuration of two parallel perfectly conducting plates, the field modes can be divided into three types. 
. Two of these are discrete modes where the potentials vanish identically in the plates, and the third type is continuum 
modes with nonvanishing potentials in the plates. The formula for the contribution of the continuum modes to the 
Casimir energy can be regarded as a special case of the Lifshitz formula 0, [H, H[ for the Casimir energy between 
■ dielectric plates. 

7-H ! I 11 @ ) we have generalized the work Q and considered the Casimir effect on a pair of magnetodielectric plates in 
an arbitrary medium. We found that the Casimir effect of a massive vector field and a massless vector field have 
significant differences. As a result, it is important to investigate the Casimir effect of massive vector fields in geometric 
configurations other than the parallel plates. In this article, we consider the Casimir effect of a massive vector field 
between two concentric perfectly conducting spherical bodies. The Casimir effect of a massless vector field on two 
perfectly conducting spherical shells have been considered in • It was shown that the Casimir force always tends 

to attract the two shells to each other. It will be interesting to see whether this is still true for a massive vector field. 



II. CASIMIR ENERGY OF TWO CONCENTRIC PERFECTLY CONDUCTING SPHERICAL BODIES 



In this section, we compute the Casimir energy between two concentric perfectly conducting spherical bodies (see 
Fig. [1} due to the vacuum fluctuations of a massive vector field. A massive vector field is represented by a four- vector 

— ,A x ,A y ,A z ), where (p is the scalar potential and A = (A x ,A y ,A z ) is the vector potential. The electric field E 
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FIG. 1: Two concentric perfectly conducting spherical bodies. 



and the magnetic field B are given respectively by 



dA 

E=-— -Vp, B = VxA. (1) 

From these, one immediately obtain two of the Maxwell's equations: 

V • B = 0, V x E+ — = 0. (2) 
at 

For a massive vector field, the other two Maxwell's equations have to be modified. In a magnetodielectric medium 
with permittivity e and permeability /i, assume that the usual linear relations D = eE and B = /iH hold. Then 

2 

„ m 

VxH -^ + 7^ A = J/ ' 

where and J/ are the free charges and the free current in the medium. Eqs. $2$ and ([3]) together constitute the 
Proca equations for a massive vector field 0] . The potentials ip and A satisfy the Lorentz condition 

due to the conservation of free charges. Both the scalar potential (f and the vector potential A have to be continuous 
across the interfaces of media @, [H, [ll[ . It has been shown in j| that a sufficient set of boundary conditions is the 
continuities of if, A and H||. 

In the following, we find the Casimir energy of a massive vector field in the configuration depicted in Fig. [TJ where 
a perfectly conducting ball of radius a\ is in the center of a perfectly conducting spherical shell occupying the region 
a 2 < r < (13. For simplicity, we only consider the scenario where the thickness of the spherical shell is infinite, i.e., 
03 — > 00. 

As discussed in 0, [B| , the eigenmodes of a massive vector field can be divided into transverse modes with V • A = 
and longitudinal modes with V x A = 0. The transverse modes can be further divided into type I modes with E r = 
and type II modes with B r = 0. They are generalizations of the TE and TM modes for a massless vector field. To 
satisfy all the boundary conditions of a massive vector field, type I transverse modes can be treated alone and they are 
called TE modes. However, the type II transverse modes and the longitudinal modes have to be considered together 
and their linear combinations are called TM modes @. 



A. TE modes 



TE modes are type I transverse modes. The continuity of H|| imply that the potentials vanish identically in the 
perfectly conducting objects. In the vacuum between the objects, the potentials are given by ip — 0, A r — and 



P, m (cos6») 
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where I = 1,2, ... , —l<m<l, 

,2 ™2„2 



9 uj m c 



A = 



ft 2 ' 

P[ n (z) are associated Legendre polynomials and ji(z) and yj(^) are the Riccati-Bessel functions defined by 

3i(z) = \Jy Ji+$(*)> m(z) = J^Y l+ i(z). 

J v {z) and Y v (z) are the Bessel functions of first kind and second kind. The continuity of A|| gives the boundary 
conditions: 

2tj i (Aa 1 ) + «B^(Aai) = 0, 
2^(Aa 2 ) + <8y ; (Aa 2 ) =0. 

Therefore, the TE eigenfrequencies are solutions of the equation 

j;(Aa 2 )yi(Aai) - ji(\a 1 )y l (\a 2 ) = 0. 

The TE modes are direct generalizations of the TE modes in the massless case [3-[To| ■ By requiring that the Casimir 
interaction energy between the spherical objects to be renormalized so that it approaches zero when the separation 
between the objects becomes infinite, one finds immediately as in [5l.[l0| that the TE contribution to the renormalized 
Casimir interaction energy is given by 



i=i 

where 



e;(7 a i) s i(7a2)' V c 2 ft 2 

Here si(z) and ei(z) are modified Riccati-Bessel functions of first kind and second kind defined by 



Iv(z) and K v (z) are the modified Bessel functions of first kind and second kind. Using the fact that t\{z) j si(z) is a 
positive decreasing function of z, we find that the TE contribution to the Casimir interaction force between the two 
bodies is always attractive. 



B. TM modes 



The TM modes are linear combinations of the type II transverse modes and the longitudinal modes. The continuity 
of H11 implies that in the perfectly conducting bodies, the type II transverse modes have to vanish. Only the 
longitudinal modes can exist in the perfectly conducting bodies. Let 



1 

c V c 2 



Inside the region r < a\ occupied by the perfectly conducting ball and the region a 2 < r < 03 occupied by the 
perfectly conducting spherical shell, the potentials of the TM modes are given by 

<p =iur- Y [Cji(Aor) +&w(Aor)] P l '(costf)^ ' l e** 



A r =- r- 2 (<£i [ji{X r) - A rj ; '(A r)] + & [yi(X a r) - A r$(A r)]) iT(«>s0)e im * 
a -2 r* h \ \i dP™ (cos 9) im x _ iuit 

P m fens 0) 

A =imr 2 [<£ji(A«,r) +&w(Aor)] ' > e^e^ 



e lut 
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Here i = 1 for the region r < a± and i — 3 for the region a 2 < r < 0,3. Moreover, 3i = 0. One can check immediately 
that the electric field and the magnetic field vanish identically in the perfectly conducting bodies. 

In the vacuum region a% < r < a 2 , the potentials are linear combinations of the type II transverse modes and the 
longitudinal modes given by 



zc 2 A 2 



[<£ 2 j/(Ar)+S 2 yKAr)]P/™(cos0)e 



A r =r- 2 {l(l + 1) [Cj,(Ar) + 2)jf,(Ar)] - <5 2 [ji(Xr) - \r~jl(\r)] - 3 2 [ Vl (Xr) - Xryi(Xr)]} If 1 (cost 
A e =r- 2 {Ar [C#(Ar) + 3)#(Ar)] + [^ 2 Ji(Xr 



\ e im<j> e -iu)t 



£ 2 y,(Ar)] } dP ™( c ° s6 ) e imt e -iu* 



=imr~ 2 {Ar [<£j,'(Ar) + Dy'^Xr)] + [£]i(Xr) + ffjr,(Ar)] } 



dO 

P; m (C0S< 

sin 6* 



(5) 



The continuities of </?, A r , Ay give rise to the following boundary conditions: 

1- ^3ji(Xoa 3 ) +S r 3 y;(A a 3 ) = 
c 2 A 2 

2. wgij'^Aoai) = [& 2 ji(Xa x ) + 

c 2 A 2 

3. u> (^3ji(X a a 2 ) +53j/;(Aoa2)) = — (£ 2 j;(Aa 2 ) + 5'2j/;(Aa 2 )) 

4- - (ii(^oai) - Aoaii;(Aoai)) £1 
= /(I + 1) [€ji(X ai ) + a)jft(Aoi)] - (ji(Aoi) - Aaij^Aai)) £ 2 - (jfi(Aai) - Aai#(Aai)) S 2 

5- ~ (j;(Aoa 2 ) - A a 2 j;(Ao<i 2 )) (£3 - (fji(X a 2 ) - A a 2 $(A a 2 )) #3 

= I(i + 1) [<tji(Xa 2 ) + S)y;(Aa 2 )] - (j;(Aa 2 ) - Aa 2 j,'(Aa 2 )) <£ 2 - (y/(Aa 2 ) - Xa 2 y[(Xa 2 )) # 2 

6. €ij;(A ai) = Aai [£#(Aai) + 2>#(Aai)] + [€aii(^»i) +3 r aw(Aoi)] 

7. € 3 j;(Aoa2) +5 , 3y;(A a 2 ) = Aa 2 [£j;(Aa 2 ) + Dy l (Xa 2 )] + [<£ 2 ji(Xa 2 ) + 3 r 2 y;(Aa 2 )] 

The eigenfrequencies of the TM modes are those u such that this system of linear equations in (£, D, £1 , <£ 2 , (S3, $ 2 , #3) 
has a nontrivial solution. Using the same techniques as in [f| fl2f . one finds that the TM contribution to the renor- 
malized Casimir interaction energy is given by 



E^ s (s) =- ]T(2Z + 1) / In AU^R, 



where 



Q l = (Qlj) is the 4 x 4 matrix 



A 1 



with 



W{ 



wi 



wi = 



Q21 



Q13 
Q23 



Q31 

Qu 



Q12 
Q22 



Q 2 i 



Q32 
Q42 

Q l 34 
Q44 



Q l (0 



/ 7aisJ(7ai) 70 1 eJ(7a i ) 
\ia 2 s[(-fa 2 ) 7a 2 eJ(7a 2 ) 



detQ'(Q 
det 



W[ W l 2 

w l 3 wi 



( 



h 2 



-si(jai 

> 

-Si(7a 2 



2 2 
m c 

2 2 



e;(7ai) 



\ 



-ej 



(7^2 )y 



ft 2 ' W *' ft 2 

Z(Z + l)s ; (7 ai)s;(7ai) /(/ + l)s ; (7 ai)e ; (7ai) N 
1(1 + l)e;(7 a 2 )s;(7a 2 ) 1(1 + l)ei(-f a a 2 )ei(-fa 2 ) 

7 2 Si(7«i)si(7oai) - 7o s i(7oai)«i(7ai) 7 2 ei(7ai)s*(7oai) - 7o s K7oai)e;(7<2i)\ 
7 2 $1 (7^2 )e; (70 02) - 7o e K7oa2)s;(7a2) 7 2 ei(7 a 2)e;(7oa2) - 7o e i(7oa2)e;(7a 2 ) J ' 
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and Q l is the 4x4 matrix 



Qo 



^Qii Q12 Q13 Qii\ 





Q31 Q32 Q33 Q34 
Qi 3 / 



(Q12Q34 Q32Q14) {Q21Q4 



43 



Q23Q41) 



Here 



/) 2 



7o = -> 
c 



s/(z) = sj(z) - ZSj(s), e/(z) = e/(z) - ze[{z). 
The expansions of det Q l (£) and det Qd(£) are given in the Appendix. 



In the massless limit m — 0, 7 = 70. Therefore, W\ = 



- To 



-ai 
a 2 



and we find that 



A TM WU1 s K7oai)eK7QQ2) 
i-i; IK J — 1 77 \~~77 \' 

e;(7oai)s; (7002) 

This is precisely the TM contribution to the Casimir interaction energy between two perfectly conducting spherical 
shells due to a massless vector field dig. 



III. DISCUSSION AND CONCLUSION 



In the case of two parallel perfectly conducting plates, the TM modes can be decomposed into discrete modes whose 
potentials vanish identically in the perfectly conducting plates and continuum modes whose potentials do not vanish 
in the perfectly conducting plates. However, for perfectly conducting spherical bodies, one cannot find TM modes 
that vanish identically inside the perfectly conducting bodies. In fact, one can check that if one impose the conditions 
that £1 — (£3 = §3 = so that the potentials vanish identically on the perfectly conducting bodies, the boundary 
conditions §5§ imply that £ = T) = £2 = $2 = 0, which means that the potentials are identically zero. This can be 
considered as a major difference between planar and nonplanar objects for a massive vector field. 
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In Fig. [2] we plot the dependence of the TM contribution to the Casimir interaction energy as a function of 02/01. 
Here we choose a% = 1cm, and we consider the cases where m = 0, m = 10~ 5 eV and m = 10~ 4 cV. In all these cases, 
we find from the graph that the TM contribution to the Casimir interaction energy is always attractive. The total 
Casimir interaction energy as a function of 02/01 is plotted in Fig. [3J Although the analytical expressions for the 
TE and TM contributions to the Casimir interaction energy are very different, TE and TM modes contribute about 
the same to the total Casimir interaction energy. This is analogous to the case of two parallel perfectly conducting 
plates. In Fig. we show the dependence of the total Casimir interaction energy on mass when 02/01 = 1.1 and 
02/01 = 1.5. It is observed that mass corrections can significantly change the magnitude of the Casimir energy. In 
fact it can be verified analytically that the Casimir interaction energy tends to zero when the mass tends to infinity. 
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FIG. 3: The total Casimir interaction energy as a function of 02/01 when a\ = 1cm. 
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FIG. 4: The total Casimir interaction energy as a function of mass when a\ = 1cm and 02/111 = 1.1 or 1.5. 



In conclusion, we have studied the Casimir effect on two perfectly conducting spherical bodies due to a massive 
vector field. In contrast to the scenario of two perfectly conducting plates, there are no discrete TM modes that vanish 
identically in the perfectly conducting spherical bodies. As a result, the analytical formula for the TM contribution to 
the Casimir energy is much more complicated. Nevertheless, numerical simulations show that the Casimir interaction 
force is attractive. 



Appendix A: Alternative expression for A^m 

The expansions of the determinants detQ(£) and detQg(£) that appear in A TM (i£) are given by 

detQ< (0 =4(0 + -jjH(* + 1)^0(0 + [-^) l 2 (l + 1) 2 CS(0, 

where 

•40 =([7 2e i(7 a i)^(7oai) - 7o s i(7oai)e;(7ai)] [7 2s i(7 a 2)e;(7oa2) - 7o e *(7oa2).?((7 a 2)] 

- [l 2 si("fai)si("/ Q ai) - 7os ; (7 ai)s z (7ai)] [7 2 e;(7a 2 )e ( (7 a 2 ) - 7o e * (7002)^(702)]) 
x ^7aie;(7ai)7a 2 s;(7a2) - 7aiS;(7ai)7a 2 e;(7a 2 )) , 

B\Q ^si(-i ai){^a 2 s' l (^a 2 )ei(^a 1 ) - "fa 2 e' l (ja 2 )s l ("fai)j 

x (ei^ai) [y 2 si (702)6; (7002) - 7o e * (7002)5; (702)] - Si(7«i) [7 2 eK702)e;(7o02) - 7o e K7oo 2 )e ( (7a 2 )] ) 
+ e;(7oa 2 )^7aisJ(7ai)e;(7a 2 ) - jaie'^a^si^a^ 

x ^(702) [7 2 s;(7ai)s;(7oai) - 7 2 s ( (7 ai)s ( (7ai)] - s ( (7a 2 ) [7 2 e;(7ai)s ( (7 ai) - 7 2 s ; (7 ai)e;(7ai)]) 

- 2j 2 j„a 1 a 2 si(j a 1 )e l (j a 2 ), 

C Z (0 =s;(7oai)e;(7 a 2 )^e;(7ai)s;(7a 2 ) - s;(7ai)e ; (7a 2 )) , 

and 

4(0 =7ai e ;(7 a i)7 a 2sK7a2)[7 2 eK7ai)S;(7oai) - 7o s *(7oai)e;(7ai)] [7 2s i(7«2)e;(7oa2) - 7o e K7oa2)s;(7«2)] , 
#o(0 =s;(7 a 2) 2 7 a i e !(7 a i) e i(7oa2)[7 2e i(7 a i)^(7oai) - 7o s *(7oai)e;(7ai)] 

+ e ( (7ai) 2 7a2sK7 a 2)s;(7oai)[7 2 s;(7a2)e ; (7oa 2 ) - 7o e i(7oa2)s;(7 a 2)] , 
c o(0 =s ; (7a 2 ) 2 e;(7ai) 2 s ( (7oai)e((7oa 2 ). 

One can show that A l (£), B l (£), C z (£), 4(0; S o(0> c o(0 are a11 positive functions of £. 
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